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Abstract 

We estimate the non-perturbative power-suppressed corrections to heavy 
flavour fragmentation and correlation functions in e^e~ annihilation, us- 
ing a model based on the analysis of one-loop Feynman graphs containing 
a massive gluon. This approach corresponds to the study of infrared 
renormalons in the large- rij limit of QCD, or to the assumption of an 
infrared-finite effective coupling at low scales. We find that the leading 
corrections to the heavy quark fragmentation function are of order A/M, 
where A is a typical hadronic scale (A ~ 0.4 GeV) and M is the heavy 
quark mass. The inclusion of higher corrections corresponds to convolu- 
tion with a universal function of M(l — x) concentrated at values of its 
argument of order A, in agreement with intuitive expectations. On the 
other hand, corrections to heavy quark correlations are very small, of the 
order of {X/QY, where Q is the centre-of-mass energy and p >2. 
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1 Introduction 



The production of hadrons containing heavy quarks in e'^e~ annihilation has proved 
to be a process of great importance for testing the Standard Model and searching 
for new physics. Heavy flavour processes are also valuable as a testing-ground for 
new techniques in QCD, because the heavy quark mass M provides a second large 
momentum scale, in addition to the overall hard scale Q, set by the centre-of-mass 
energy in e+e~ annihilation. This has led to a rather good understanding of heavy 
quark production and fragmentation within the context of perturbation theory [0-^]. 
On the other hand, it has become clear as a result of this understanding that the c and 
b quark masses are not large enough for a purely perturbative approach to provide a 
good description of the data for these flavours [||,|| . This is because non-perturbative 
effects can give rise to contributions of order A/M, where A is a typical soft physics 
scale. One may also worry about the possibility of non-perturbative effects of order 
X/Q, which could be significant at present energies. Thus it becomes important to 
study power-suppressed corrections in QCD from as many viewpoints as possible, and 
to apply the resulting insight to heavy flavour processes in particular. 

One approach to the study of power-suppressed corrections, which has proved pop- 
ular recently, arose from the study of infrared renormalons . Here a divergent series 
of perturbative contributions gives rise to a power-suppressed renormalon ambiguity 
in the prediction of perturbation theory. One can then argue that a non-perturbative 
contribution with the same power behaviour should be present, with an ambigu- 
ity in its coefficient which cancels that associated with the renormalon. Conversely, 
when the renormalon ambiguity involves a high inverse power of the hard scale, one 
expects non-perturbative corrections to be especially small. This approach can be 
reformulated without reference to renormalons by postulating the existence of an 
infrared- regular effective coupling at low scales (the 'dispersive approach' of ref. 0). 
The expected power-suppressed corrections are consistent with those predicted by 
more rigorous approaches such as the operator product expansion where applicable, 
and have been found to agree fairly well with those suggested by experimental data 

The aim of the present paper is to study power suppressed effects in heavy flavour 
production in e~^e~ annihilation, using either the 'renormalon' or 'dispersive' ap- 
proach, which are equivalent for the purpose of the present work. Our aim will be to 
test some standard assumption about the form of non-perturbative corrections that 
enter the fragmentation function. Furthermore, we will also examine the correlation 
between the quark and antiquark momenta. This topic is of practical interest, since 
this correlation affects the determination of r^, the partial width for the decay of the 

boson into 6-flavoured hadrons. 
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1.1 Renormalons and power-suppressed corrections 

In the renormalon approach, described in detail for e^e~ shape variable calculations 



in ref. |]T3[, one starts from the first-order perturbative contributions to the process 
e^e~ QQX, involving emission of a single real or virtual gluon. One then identifies 
a factorially divergent series of perturbative contributions associated with light quark 
loop insertions on the gluon line. These contributions will be dominant at high orders 
for sufficiently large values of n/, the number of light flavours. One now argues jl^ 
|rB[ that the true high-order behaviour of perturbative QCD can be approximated 
by making the replacement Uf ^ nf — 33/2 in the large-nj behaviour. There is 
support for this 'naive non-Abelianization' assumption in the high-order behaviour 
of the average plaquette in quenched {nj = 0) lattice QCD [ p!7[] . 

The resulting factorial divergences of the perturbation series can be of two types. 
Those associated with high-momentum regions of integration (ultraviolet renormalons) 
correspond to divergent series with alternating signs, which can be summed unam- 
biguously using standard techniques such as Borel summation. We shall not be con- 
cerned with them in the present paper. Those due to low momenta flowing in loop 
integrals (infrared renormalons) produce same-sign asymptotic series, which are in- 
trinsically ambiguous. The ambiguity is of the order of the smallest term in the 
series, which turns out to be a power-suppressed quantity, of order {\/QY where p is 
a number that depends on the observable being computed. 

In the full theory of QCD, the infrared renormalon ambiguities of perturbation 
theory must be cancelled by non-perturbative contributions. We shall assume that 
the presence of an infrared renormalon with a particular value of the power p indicates 
that a comparable power-suppressed non-perturbative contribution is actually present 
in the full theory. 

The dispersive approach of ref. involves similar calculations, with a slightly 
different interpretation. A (formal) dispersion relation for the QCD running coupling 
of the form 

«.(^') = - r ' = -IT-- Di^^ {"^( V)} , (1-1) 

is assumed. In QCD, unlike QED, the running of the coupling cannot be associated 
with vacuum polarization effects alone. However, in the same spirit as the 'naive non- 
Abelianization' assumption, it is further assumed that the dominant effect on some 
QCD observable F of the running of as in one-loop graphs may be represented in 
terms of the spectral function PsifJ^"^) and a characteristic function J^{ij?), as follows: 



a,(0)^(0)+ r^p,(^2)._^(^2 

Jo L 
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where the relation ( |1 . 1| ) has been used to ehminate as{0). The characteristic function 
is obtained by computing the relevant one-loop graphs with a non-zero gluon 
mass fi [p!^p| and dividing by as- Integrating Eq. ( p,.2D by parts, we can write 



/i^ 



aefr(/i') ■ ^(/i') , (1.3) 



where 



and we have introduced the effective coupling aesifi"^), defined in terms of Ps{l^'^) by 

It follows from this definition and Eq. ( |1 . 1| ) that 

ae.(/. ) = «.(/. ) - + -... = «. + 0(aj , (1.6) 

and therefore in the perturbative domain <^ 1 the standard and effective cou- 
plings are approximately the same. In the large-n/ renormalon approach^] we have 
the explicit expression, obtained by substituting the one-loop running coupling into 
Eq. (^: 

aofr(yU^) = — ;— arctan ( 7r6oas(/i^) ) (1-7) 

where 6o = (33 — 2n/)/127r. 

The characteristic function Ti^^j?^ is more precisely a function of the dimensionless 
ratio e = p'^ /Q"^, where is the characteristic scale of the hard process. If the effec- 
tive coupling in Eq. (|l]3|) has a non-perturbative component 5Q;efj(/i^), with support 
limited to low values of /i^, the corresponding correction to F, 

5F= f ^Sa,^{^^') J'ifi') , (1.8) 

will therefore have a dependence determined by the low-/i^ behaviour of JF. 

A crucial point is that only those terms in JF that are non-analytic at /i^ = (e = 0) 
can produce power-suppressed contributions to 6F. This is because the integer /z^- 
moments of (5aefr(/^^) are required to vanish, for consistency with the operator product 



expansion. The same result may be seen in the renormalon analysis of ref. [13|: for any 



behaviour of JF(e) of the form e*^ as e — > 0, the renormalon contribution is proportional 

^We stress that this model is not physically fully consistent, because of the presence of the Landau 
pole, which implies that the support of the spectral function of must be extended to negative 



values of the argument for Eq. (1.1) to be valid 
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to Im(e*^'^) and therefore vanishes for integer p. On the other hand, a -y/e behaviour 
imphes a non- vanishing correction proportional to 1/Q, while elne gives etc. 
Thus our objective is to identify the leading non-analytic terms in the behaviour of 
the characteristic function at small values of the gluon mass-squared, which will tell 
us the Q^- dependence of the leading power-suppressed corrections. 

If one makes the additional assumption that the effective coupling modification 
(5aefr(/i^) in Eq. ( p..8D is universal, one obtains a factorization property for power- 
suppressed corrections, which leads to relationships between the coefficients of the 
corrections to different observables. For variables like event shapes, this type of 
factorization is only approximate, due to the fact that a cut dressed gluon line is 
weighted differently, depending upon the value of the shape variable for the particular 
final-state structure of the cut gluon [|1^]. In the present case, however, the dressed 
gluon is cut fully inclusively, without any weight, and therefore factorization may be 
more reliable. 

In the case of heavy fiavour processes, we also want to study corrections that are 
suppressed by powers of the heavy quark mass, M. As long as we treat both M and Q 
as large parameters, and keep track of the dependence on their ratio, this will be done 
automatically when we extract the non-analytic terms in e. Defining p = AM'^/Q'^, a 

y^e/p term will indicate a correction proportional to 1/M, elne/p implies and 
so on. 

In our terminology, mass corrections of the form {M/QY will not be called power- 
suppressed, since we are always assuming that M is not small. 



2 Calculations 



2.1 Massive gluon cross sections 

Considering first the vector current contribution, the distribution of the heavy quark 
and antiquark energy fractions x and x with emission of a gluon of mass n in the 
process e'^e~ QQg is given by 



1 (fay as Cp 



ay dx dx 271 (3 



V 



V 



{x + rif + (x + rif + Cv 

(l + ip)(l-x)(l-S) {l-xf {l-xf 



(2.1) 



Here 1] = {fi'^ + 2M^)/Q^ = e + |p where e = p^/Q^, M is the quark mass, p = 

Cv = -2p(l + T]) , (2.2) 



P = a/1 — p is the heavy quark velocity, and 

av = ao{l + ip) P (2.3) 

is the Born cross section for heavy quark production by a vector current, ctq being 
the massless quark Born cross section. 
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The phase space is determined by the triangle relation 



A{x'^ - p,x^ - p,xl-4e) <0 (2.4) 

where A(a, b, c) = + 6^ + — 2ab — 2bc — 2ca and Xg = 2 — x — x. This gives 
x^ < X < x^ where 

_ (2 - x)(l - X - e + |p) ± p, e) 

^ 2il-x) + lp (2.5) 



S(x,p,e) = ^(x2-p)[(l-a;-e)2-ep] , 

and 

v/p<a;<l-e-yip. (2.6) 
In the case of the axial current contribution, instead of Eq. ( |2.1| ) we have 

1 (PaA _ as Cp \{x + r]Y + {x + + Ca V V 

a A dx dx 271 [3 

where 



(l-p)(l-a;)(l-x) (l-x)2 {l-xf 



(2.7) 



a = ip[(3 + a;,)2-19 + p-8e], (2.8) 
a A being the Born cross section for heavy quark production by the axial current: 

aA = ao/?' = ao(l-p)/5. (2.9) 
2.2 Leading power corrections 

Clearly the expressions ( |2.1j ) and ( |2.7| ) are analytic functions of e at e = except 
possibly for x and x near 1. The phase space is also analytic in e whenever the gluon 
momentum is large, since in this region one can always expand kinematic variables 
in powers of e. As discussed above, this implies that there are no non-perturbative 
corrections of the type we are considering for x, x < 1. Non- analytic behaviour may 
only arise in the region where Xg ~ (x, x ^1). In order to investigate the 
corrections associated with this region we take moments of the form 

M{N,N,e) = I dxdxx^'x^--^ (2.10) 
J a dxdx 

and expand 

x^x^ = 1- Ny- Ny + --- , (2.11) 

where y = 1 — x and y = 1 — x. The first term corresponds to the total heavy 
flavour cross section, whose dominant power correction is of order {\/QY or smaller. 
However, the next two terms give contributions proportional to y/e at small e, which 
could give rise to \/Q and/or A/M corrections. To evaluate them we note that their 
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contribution to the difference 5Ai = Ai {N, N,0) — (N, N, e) for small e can be 
written, for both the vector and axial current contributions, as CpO-s^ /'^'^ where 



T = ^ Jjydy{Ny + Ny) 



V yy 2y 

'4 - 3p py 



13 J R \ yy 2y2 2y 



r,2 



1-(N + N) f dydy[ 
213 Jr \ y y 



(2.12) 



R being the region between the phase space boundaries for e = and e > 0. Eq. ( |2.4| ) 
may be expanded in the region Xg = y + y ^ y/e, so that the boundary of phase space 
in this region is given by 

iA(x2 ~p,x'- p, xl - 4e) ^ p{y^ + f) - 2(2 - p)yy + 4e(l - p) < , (2.13) 

which is the equation of a hyperbola. Changing variables to r, where y = r cos 4> 
and y = r sin 0, the region R may be written as 



R: S<4.<h'-S, 0<r<J^p^^ (2.14) 
^ V sm 2(p — sm 2d 

where sin 25 = p/{2 — p). Performing the integration one gets 



^ 71" /2e(l - sin2(5) „ p 
= —{N + N)d ^ 4-3p ^ 



2(3^ 1 sin 25 V sin 25 

2.^^(iV + iV)7I=.(iV + iV)(l-4^)A 



;2.i5) 



Therefore the leading power correction is of order p/M rather than p/Q. The leading 
correction with an explicit dependence on Q is of order pM/Q"^. This is consistent 



with the finding for light quark fragmentation functions ||g,|T^,0,0: the leading 
power correction in Q is of order l/Q"^. 

The linear dependence on the moment index in the result ( p.l5| ) implies a 
behaviour in x-space of the form 

J^{x, X, p^) = 7Tp^^ [5(1 - x)5'(l - x) + 5'(1 - a;)5(l - x)] . (2.16) 

This means that, as far as the leading power correction is concerned, the two-particle 
heavy quark distribution factorizes. In the dispersive approach of ref. 0, the non- 
perturbative correction is given in terms of the low-energy modification to the effective 
coupling, Saeffip'^), by Eq. ( |1.8| ). Defining the non-perturbative parameter 



Ai = ^ j^p6a,s{p'), (2.17) 
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Eqs. ( |l.(j| ) and ( 2.16 ) then imply that 



1 £a 
a dx dx 



6{1 



x] 



TcAi 

2M 



p^6\l - x) 



5[l-x-p' 



M 



■6 1 




2M 



(3'S'{1 - x) 



(2.18) 



where A = ttAi/2. Thus we see that the main non-perturbative effect is a shift in the 
heavy-quark momentum fractions by an amount 6x ~ X/M. Assuming approximate 
universality of Sacs, one may estimate from hght-quark event shape data that Ai ^ 
0.25 GeV [0, which gives A ~ 0.4 GeV. This agrees with the order of magnitude of 
the non-perturbative shift estimated from (x) in heavy flavour fragmentation [0. 



2.3 Higher power corrections 

Power-suppressed effects in the heavy flavour fragmentation functions should be 
equivalent to a convolution with a non-perturbative initial condition of the form 
Mf (M(l — x)), and therefore should approach a delta function as M — >■ oo. This 
can be inferred by intuitive reasoning, but can also be derived more rigorously in the 



context of the heavy quark mass expansion |jT9[. In this section we will show that 
this expectation is also fulfilled in our model. 

First of all we note that Eq. (|2.18|) can be written as 

1 d^a 



a dx dx 



M6 (m(1 - x) - P^X) ■ M6 (m(1 - x) - P^X) (2.19) 



and therefore the expected form is indeed obtained when one includes only the leading 
power correction. 

To go beyond the leading correction, we have to consider higher moments with 
respect to y and y in Eq. (|2.12|) . We examine first the moments of the single-particle 



distribution (A^ = 0). For moments weighted by y^ with p > 1 we have to define 
the integration region R more carefully. Consider for simplicity the case that p is 
small (i.e. -C Q"^), so that sin 25 ^ 26 ^ |p. The upper limit of the r-integration 
becomes a constant of order unity when S < (p < 6 + e. Hence this region gives a term 
that is analytic in e, which will not contribute to power corrections. The important 
region is 6+e < (p < |7r— 5 — e. The leading non-analytic term coming from this region 
is proportional to (e/5)^/^ when p is odd, and proportional to (e/5)^/^lne when p is 
even. Hence for every value of p there is a power correction of order (e/(5)^/^ = {n/My. 
In detail, for p small and iV = Eq. ( |2.12| ) becomes 




R '''' \yy 2?/2 



dydyy 

rl 

=+5 0' 



(2.20) 

Je+5 0'^ 
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Now 

^(0 _ <5)(2-P)/2 ^ (_l)(p-l)/2 _ O ^5-p/2 ^ . . . 

02 ^ V 2 y ^2.21) 

~ (_l)p/2 _ Ip^ 5-P/2 ine + . . . for p even, 

where the dots correspond to terms giving contributions that are analytic and/or 
higher-order in e. Thus, keeping only the leading non-analytic parts, we find 

which corresponds in x-space to 

^(x,/i2) = M/ {M{l-x)) 

f{z) = Re I (In 6 - zvr) E ^ " l) (^/^)' ^ (^) ^ ^^'^^^ 

We can therefore see the expected scaling of the fragmentation function. We also see 
that while the leading power correction, eq. ( |2.18| ), corresponds to a simple shift in 
the value of x, this behaviour is not preserved by the higher power corrections. 



3 Correlations 

We consider now the higher power corrections to the two-particle distribution, i.e. 
the inclusion of higher powers of /i/M in Eq. ( |2.1(j| ). We have to examine the double 
moments corresponding to Eq. (|2.12|) with general weights y'^y'^. Again we treat only 
the case of small p. Then we find that for any q such that < g < the leading 
non-analytic term is suppressed by a factor of p'^ relative to that for weight y'^^'^. 
Therefore the leading terms in each order of p/M remain of the form (|2.16|) : 



Tix, X, /i^) = 5(1 — x)J-'{x, p^) + J-'{x, p^)5{l — x) (3.1) 

where T{x,p^) is as given in Eq. ( |2.23| ). Thus the two-particle distribution including 
these terms still factorizes and can be expressed as a function of M(l — x) and 
M(l — x), although beyond leading order in 1/M it differs from the simple product 
of delta functions given in Eq. ( p. 19 ). 



Because of its possible impact on the determination of Ff, in decays, it is inter- 
esting to determine what is the leading power correction to the momentum correlation 
{yy) — (?/)2. As stated before, the correction to {yy) behaves as p{^JT/pY = e at small 
e. The term (y)"^ gives zero at the order we are considering. In fact in the large- 
Uf limit it is subleading, and in the dispersive approach it is of second order in the 
effective coupling. Therefore the momentum correlation is of order e. 
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In order to confirm this conclusion, we also calculated the difference {yy) — {yy)^=Q 
using the exact phase space and matrix elements. We found that the leading term at 
small e is proportional to e In p. Thus, corrections to the momentum correlation in our 
model are suppressed by at least two powers of Q, and should therefore be completely 
negligible at LEP energies. Whether this result survives higher-order corrections is an 
open question, and in fact a very difficult one. We simply point out here that, while 
Monte Carlo models seem to indicate the presence oi \/Q corrections to correlations 
(see ref. ^ and references therein), the simple model that we have adopted in this 
work does not provide support for the presence of such corrections. This is also 
consistent with the findings of ref. [|1^], where power corrections to fragmentation 
functions were computed in the strictly massless limit. 



4 Discussion 

We have examined the heavy flavour fragmentation function and correlations in a 
simple model, and found the following results. 

At leading order, non-perturbative effects in the fragmentation function can be 
represented as a convolution with a function of the form 

M/(M(l-a;)), (4.1) 

which approaches a 5 function as M oo. The leading power correction has the form 
A/M, where A is a typical soft hadronic scale. An estimate of this correction based 
on the approach proposed in ref. |0 gives the correct sign and order of magnitude. 

In the two-particle distribution, corrections of the order of (A/M)^ factorize and 
therefore no large correlations of this order arise. Since correlations are important 
for their possible impact on the determination of Ti, in decays, and since the 
perturbative value of the correlation is of the order of 1% [Q, it is also important 
to understand whether corrections of the order oi \/Q are present. In our analysis, 
consistently with ref. [0, terms of this order do not arise. In fact, we also verified 



numerically that the e dependence of the correlation {yy) — (y)"^ is of order elogp, 
and therefore the leading power correction is less than order (A/(5)^logp. 

We end with a comment on the relationship between our results and those of 
ref. 0. In that paper a resummed expression for the heavy quark spectrum was 
derived and numerical results were presented using various models for the behaviour 
of the QCD running coupling at low scales. The non-perturbative component of the 
coupling generates 1/M corrections which should correspond to those considered here, 
after convolution with the perturbative fragmentation function. 
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Appendix 

We give here for reference the single-particle inclusive distribution which results from 
integrating the vector two-particle distribution (|2.1|) for x < 1: 



1 day dx d^ay 



CTy dx Jx- cry dx dx 



X In 



27r/5(l + |p)\V l-x 

[(1 -x){x- ip) + e(2 -x) + p, e)f 



(A.l) 



(1 - X + ip) [p(l - x)2 + 4e(x + e - p)] 

S(x,p,e) $(x,p,e) 

]{1 - - x + ip)2[p(l - x)2 + 4e(x + e - p)] 



where a;± and S are given in Eq. (2^) and 



$(x,p,e) = p(l-a;)2[4(l-x)2(8-x) + 2p(l-x)(17-92;) + 2p2(5-4a;)+p3] 
+ 2e[8(l - xf{2 + x'^) - 2p(l - x)(4 - 7x - 12x2 + 7x^) 
- 2p2(13 - 19x + + 4x^) - p^(9 - 7x - x^) - p^] 
+ 2e2[8(l - x)(4 - X - 2x2) + 4p(l + 5x - Sx^) - 10p2(l - x) - p^] 
+ 4e^[4(l - x)(4 - 3x) + 2p(5 - 4x) + p^] . 

(A.2) 
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